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Low-energy proton captnre reactions in the mass region 55-60 
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Low energy proton capture reactions in the mass 55-60 region are studied in a microscopic optical 
model. Nuclear density profile is calculated using the relativistic mean field theory. The DDM3Y 
interaction is folded with the theoretical density to obtain the proton-nucleus optical potential. 
A definite set of normalization parameters has been obtained for the concerned mass region by 
comparing with all available experimental data in this mass region. These parameters have been 
used to obtain proton capture rates for astrophysically important reactions in this mass region. 

PACS numbers: 24.10.Ht, 25.40.Cm, 25.40.Lw, 25.40.Kv 


I. INTRODUCTION 

The study of the role of individual nuclear reactions 
in stellar evolution has been an important field of re¬ 
search in the last few decades. As a star evolves with 
time it passes through burning in different ranges of nu¬ 
clear mass. At the same time, different nuclear processes 
become important at different time periods of evolution. 
A comprehensive study of these processes sheds light on 
various astrophysical phenomena. 

There are certain astrophysical sites which are respon¬ 
sible for the production of heavier nuclei beyond iron 
through the rapid capture of protons on seed nuclides. 
In the mass region of our interest there are certain pro¬ 
ton rich naturally occurring nuclei, which are not pro¬ 
duced by the r-process or the s-process. These are called 
p-nuclei. Proton capture reactions in certain astrophys¬ 
ical sites can account for the formation of some of these 
proton rich nuclides. For example x-ray bursters with 
a large proton flux in the peak temperature around 1-3 
GK are suitable astrophysical sites for the production of 
certain nuclei. To find out the abundance of different 
nuclei as well as the evolution of the process in these 
sites a network calculation is necessary which involves a 
large number of reactions. 

It is thus imperative to calculate the rates and/or cross 
sections of these reactions in different mass ranges. Our 
group has already calculated the cross sections and hence 
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the astrophysical S-factors in the mass range A = 60 — 
100 [IH|. Some implications of the new rates has also 
been investigated in the context of rp-process [ 1 , @ . In 
the present work, we extend our calculation to the A = 
55 — 60 region. 

The rp-process is sensitive to a number of reactions 
in this region. The most challenging aspect to look at 
in these scenarios is that most of the nuclei involved in 
those reactions are not produced in the laboratory. For 
example, Parikh et al. [7| have identified proton capture 
reactions on ®®Ni and targets as important in the 

rp-process in certain scenarios. However, experimental 
rates are not available for these reactions because stable 
targets do not occur in nature. Hence, one has to depend 
on theoretical calculations in this domain. 

In explosive proton rich environments, such as x-ray 
bursts, proton capture has to compete with its inverse, 
i.e. photo-disintegration. This competition results in 
waiting points and causes delay of further nucleosynthe¬ 
sis. With temperature, pressure and proton mass frac¬ 
tions being different at different regions of these sites as 
well as being time-varying quantities, incorporation of 
all these physical conditions in the nuclear network is a 
big challenge. Rauscher et al. Si have calculated the 
rates for various proton, neutron and a-particle induced 
reactions and their reverse reactions in Hauser-Feshbach 
formalism for targets with wide range of atomic numbers 
and masses and for a wide range of temperature. Theo¬ 
retical calculations in this mass region essentially ntilize 
the Hauser-Feshbach formalism where, the optical model 
potential, a key ingredient, is often taken in a local or a 
global form. However, a more microscopic approach is 
also possible using an optical potential constructed uti¬ 
lizing nuclear densities. If the target is stable, nnclear 
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density is available through electron scattering. How¬ 
ever, in the absence of a stable target, theory remains 
our sole guide to describing the density. It is impera¬ 
tive to test the theoretical calculations, where experi¬ 
mental data are available, to verify its applicability. We 
aim to check the success of microscopic optical potentials 
based on mean-field densities in explaining the available 
reaction cross sections in this mass region. A good de¬ 
scription depending essentially on theory will allow one 
to extend the present method to the critical reactions, 
which are beyond present day laboratory capabilities. 

A well defined nucleon-nucleon (NN) interaction is of 
major importance for microscopic calculation of nucleon- 
nucleus and nucleus-nucleus potentials used in the the¬ 
oretical analysis of different reactions as well as scat¬ 
tering. The optical model potential is highly successful 
for explanation of different branches of nuclear reaction. 
It can reliably predict the basic observables such as to¬ 
tal and partial cross sections, elastic scattering angular 
distributions, etc, even for those target nuclei and for 
those energy regions for which no experimental data ex¬ 
ist. We have used the density dependent M3Y interac¬ 
tion by folding the potential with target radial matter 
densities. This interaction has been used in many calcu¬ 
lations and has given satisfactory results. 

The paper is organized as follows. In the next section, 
we outline our method of calculation. Essentially we 
construct an optical model potential through folding an 
NN interaction with the theoretical density profile. For 
this purpose we use the relativistic mean field (RMF) 
theory to obtain the density profile of the targets. In 
Sec.lm] the results of our work are discussed in detail. 
Finally we summarize our work. 


II. MODEL CALCULATION 

The RMF approach has proved to be very successful in 
describing various nuclear properties such as binding en¬ 
ergy of nuclei in ground states as well as excited states, 
nuclear density profile, rms charge radii, deformation, 
nuclear halo, moment of inertia, etc M- It is considered 
to be the relativistic generalization of the non-relativistic 
models such as Gogny force or Skyrme force Hartree- 
Fock theory using effective mesonic degrees of freedom 
rather than instantaneous forces. The model is basically 
based upon two major approximations namely mean- 
field approximation and no-sea approximation [llj . The 
starting point of RMF is a suitable Lagrangian density 
that includes the coupling between the nucleon field and 


meson field as well as meson self couplings so that the 
Lagrangian can successfully describe the properties of 
finite nuclei as well as the equation of state (EOS) of 
nuclear matter. There are different variations of La¬ 
grangian density as well as different parameterizations. 
An accurately calibrated relativistic Lagrangian density, 
FSUGold [I 4 I, has been fitted to the charge radii of nu¬ 
clei. It contains two additional parameters, compared 
to conventional RMF models, describing self coupling of 
vector-isoscalar meson and coupling between the vector- 
isovector meson and vector-isoscalar meson. These two 
additional parameters significantly affect the softening of 
the EOS, the accurate determination of which is needed 
for the study of various nuclear properties such as charge 
radii, masses, etc. 

Thus theoretical density profiles are extracted in the 
RMF approach considering the FSUGold interaction. 
The charge density is obtained by convoluting the point 
proton density considering the finite size of the nucleus. 

Pch{r) = e j p(r/)g(r - r/)dr/ (1) 

where g{r) is the Gaussian form factor given by, 

gir) = {ay/^)~^exp{-r‘^/a?) (2) 

where a is a constant whose value is assigned to 0.8 fm. 
Using the nuclear density profile we have numerically 
obtained rms charge radii. 

While calculating the charge density or the radius, no 
attempt has been made to take the correction due to cen¬ 
ter of mass into account. Galculations on harmonic os¬ 
cillator wave functions show that the correction is small 
for heavier nuclei. For example, Quentin has shown [T3l| 
that the effect of inclusion of center of mass correction 
in the radius is given by 5r/r ~ 0.9/A^/^. Hence, we do 
not expect the density profile to be affected significantly 
due to this approximation. 

The M3Y interaction [l3, [il| is based on a realistic 
G-matrix which in turn is constructed in a harmonic os¬ 
cillator representation averaging over a range of energies 
as well as densities. It has no explicit density depen¬ 
dence nor energy dependence. Although in most cases 
these averages do not matter producing satisfactory re¬ 
sults, in few cases it becomes necessary to incorporate 
explicit density dependence into M3Y interaction and 
then it is named as density dependent M3Y (DDM3Y) 
effective interaction [l^ . Low energy proton capture re¬ 
actions are highly sensitive to nuclear radius as well as 
density. In the present work we have used density de- 
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pendent M3Y Reid-Elliot effective nucleon-nucleon in¬ 
teraction within a folding model prescription 0. The 
density dependence is incorporated in the same way as 
suggested in Refs. [13 [ 3 . 

Further, we have included a spin-orbit term into 
the potential considering Scheerbaum prescription 3 
which has been coupled with the phenomenological com¬ 
plex potential depths. These depths are functions of en¬ 
ergy which are assigned standard values as in Lahiri et 
al. [ 3 . These values are kept unaltered throughout our 
present work. 

We have incorporated the density dependent M3Y in¬ 
teraction within the TALYS1.4 code [3 performed 
a Hauser-Feshbach (HF) calculation. We have chosen 
Goriely’s microscopic level densities and Hartree Fock 
Bogolyubov model for El 7 -ray strength function. As 
seen in our previous calculations Si , these choices can 
explain the experimental results more accurately. All 
these options are available in the code. We have also 
included the effect of the width fluctuation correction 
which has a significant impact at low incident energies. 
Up to 30 discrete levels are included for both target and 
residual nuclei, which are considered in Hauser-Feshbach 
decay and 7 -ray cascade. We also include a maximum 
of 30 discrete levels for the nuclei resulting from binary 
emission in Hauser-Feshbach decay and 7 -ray cascade. 
HF calculations are done with full j, I coupling. We have 
incorporated the density data obtained from RMF ap¬ 
proach to obtain the optical model potential. 

Because of rapid variation of cross-section with energy 
in the low energy region, it is difficult to compare the 
theory and experiment. A standard alternative way is 
to compare another important quantity instead of cross- 
section, namely astrophysical S-factor (^|. 

The proton capture reactions in astrophysical sites oc¬ 
cur within a narrow energy window 21|. This effective 
energy window approximately of Gaussian shape around 
a peak (known as Gamow peak) is known as Gamow win¬ 
dow. The expressions for the Gamow peak and Gamow 
width in a practical form are given respectively as, 


Eo = 0.1220{Z^Zl^T^Y/^MeV 


(3) 


A = 0.2368{Z^Z^^lTiy/^ (4) 

where! /r is the reduced mass and Tg denotes the temper¬ 
ature in GK. Thus most of the astrophysically important 
reactions occur within a narrow energy window Fg — A/2 
to Eq + A/2. We see that for (p, 7 ) reactions on stable 
isotopes in the mass range 55-60, the Gamow window lies 


between 1 MeV to 3 MeV for temperature around 3 GK. 
Hence, we have carried out our calculation in this low 
energy window and compared our results with the mea¬ 
sured data where available. In calculating the Gamow 
peak, Gamow width and hence Gamow window we have 
taken the masses from Audi et al. [^ . 


III. RESULTS 

Because the optical model is dependent on the density 
profile of the nucleus, we calculate the density and the 
charge radii of nuclei in this mass region using RMF for¬ 
malism. The theoretical density values are plotted as a 
function of radius and compared with available experi¬ 
mental values in Fig. [TJ As can be seen the agreement is 
extremely good. The experimental data are taken from 
Wohlfahrt et al [^ . 

We also compare all the available rms charge radii 
values with theoretical results in table ID The exper¬ 
imental values are taken from Angeli et al. [^. It 
can be seen that the RMF calculation has an excel¬ 
lent predictive power, the relative difference between 
theory and experiment in all cases being less than 0.5%. 


TABLE I: Charge radii of various nuclei extracted in the RMF 
approach compared with measured values from Angeli et al. |24| 


Nucleus 

Charge radius(fm) 
Theory Experiment. 


3.7057 

3.7057 

®®Mn 

3.7189 

3.7146 

®®Fe 

3.7361 

3.7377 

®"Fe 

3.7497 

3.7532 

CO 

3.7634 

3.7745 


3.7924 

3.7875 


3.7916 

3.7757 

®°Ni 

3.8193 

3.8118 


We have tried to set a definite normalization for the 
optical model potential that fits all the reaction data in 
the concerned mass region. The potential obtained by 
folding has been multiplied by the normalization con¬ 
stant 2.0 to get the real part of the potential. The 
DDM3Y interaction does not have any imaginary part. 
We have multiplied the folding potential by the normal¬ 
ization constant 1.4 to obtain the imaginary part of the 
optical potential. These final parameters have been ob- 
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FIG. 1: Theoretical density profile of various nuclei in mass 55-60 region from RMF approach compared with experimental data obtained 
by elastic electron scattering taken from Wohlfahrt et. al. 12.111 . The solid lines denote the theoretical result and discrete points with 
error-bars represent the experimental data. In most of the cases, especially for lower radii values, the errors associated are smaller than 
the dimension of the empty circles. 



FIG. 2: Comparison of theoretical astrophysical S-factors calculations in the present work with experimental data for proton capture 
on (aJffMn, (b)®®Fe, (c)®®Co, and (d)®®Ni, (e)®''Ni. See text for the explanation of different symbols and the Supplemental Material in 
Ref. [2a for numerical values of S-factor. 


tained after many trials to ensure a reasonable agreement 
with experimental data for all the known low energy pro¬ 
ton capture reactions in the mass region of our interest. 
Although a single normalization can not reproduce the 


experimental data excellently for all reactions in the re¬ 
gion, i.e. each individual reaction may have different 
normalization for best matching with measurement, it is 
necessary to consider a single definite normalization to 
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extend the work to unknown nuclei in the mass region 
for which no experimental data exist. We note that the 
fitted parameters for the present mass range differ from 
the neighbouring mass region in our earlier calculation. 
This is possibly due to the fact that the mass selected 
in the present calculation is lighter than our previous 
regions. Possibly, the larger depths of the potential is 
required to adjust for the low mass region. 

The comparison of S-factors obtained after incorpo¬ 
rating DDM3Y interaction using the above normaliza¬ 
tion constants with experimental data are shown in fig. 
[2j The numerical values of the S-factors and the re¬ 
action rates are given in Supplemental Material [2^. 
The experimental data are taken from Ref [26l - E^ for 
^®Mn, ®®Fe and ®®Co, respectively. For ®®Ni and ®°Ni, 
experimental data are taken from Refs. [29l - l3Tj | and 
Refs. [13, m, 113, respectively. In many cases the ex¬ 
perimental data are very old. Errors are also not avail¬ 
able in some cases. For ®®Mn, ®®Fe and ®®Co, circles 
represent the experimental data. For ®®Ni, trian gles , 
squares and circles represent the data from Refs. [29l - [^ . 
respectively. For ®°Ni there are three different sources of 
data [30, HO, HO] which are denoted by squares, circles 
and diamonds, respectively. 


In all cases the solid line denotes the theoretical 
DDM3Y result. In ®®Mn, there are certain ambigui¬ 
ties in experimental data, especially in the energy range 
between 1.3 to 1.6 MeV. The experiment was done us¬ 
ing Ge(Li) detector by integrated beam current method 
more than three decades ago. However, errors are not 
associated with most of the data points. Only four data 
points in the energy range of our interest have errors 
associated with those. 

Our calculations give an excellent description of ex¬ 
perimental data for ®®Fe. The experiment for ®®Fe was 
done using Ge(Li) detector and the data was compared 
with statistical model predictions 27|. For ®®Go again 
there are large fluctuations in experimental data. Butler 
et al. [13 stated that they had observed several reso¬ 
nances in the reaction ®®Go(p, 7 )®°Ni but the resonances 
were too close to be resolved clearly. 

Our calculation for ®®Ni overpredicts the measure¬ 
ment of Tingwell et al. 3l| by a factor of ^ 2.5, ap¬ 
proximately. This experiment was carried out by both 
beam current integrated method and single target irra¬ 
diation method using Ge(Li) detector. Tingwell et al. 
also compared their data with statistical model calcula¬ 
tions. They found that their statistical calculation over¬ 
estimates the measurement by a factor of ^ 2.5 for ®®Ni, 
which agrees with our results. Cheng et al. [13 also 


measured cross section for this reaction using the acti¬ 
vation technique. Except in the energy range ^ 1.4-1 .8 
MeV, where the measurement itself has large discrepan¬ 
cies, the data agree more or less well with our theoretical 
calculations. 

For the reaction ™Ni(p, yj^^Cu, Tingwell et al. them¬ 
selves compared the experimental results with the sta¬ 
tistical model predictions and showed that normalizing 
the optical model imaginary well depth for Ni isotopes 
by a factor of 1.5 leads to a better agreement between 
theory and experiment [33l |. Our calculation, in the case 
of ®°Ni, overpredicts the experimental data of Tingwell 
et al. [^ by a factor ~ 1.5, whereas it underpredicts the 
data of Krivonosov et al. [^ by a factor ^ 0.35. 

With the above normalization, we have calculated 
the rates for 7 ) reactions identified as important by 
Parikh et al. Q • The calculated rates are compared with 
NON-SMOKER [I, i, Ull rates. The NON-SMOKER 
results are from a HF calculation based on masses from 
experimental measurements and calculation in the Fi¬ 
nite Range Droplet Model [13. Other details of the 
calculation can be obtained from the references. The 
results have been plotted in figure 3. We see that in the 
range 1-4 GK, the NON-SMOKER results differ from the 
present calculations significantly. Eor the ®®Ni(p, 7 ) re¬ 
action, although the results agree at low temperature, at 
higher temperature the NON-SMOKER rates are larger. 
On the other hand, for the other two reactions, viz. 
®^Cu(p, 7 ) and ®®Cu(p, 7 ), our calculation predicts a sig¬ 
nificantly larger rate throughout the temperature range. 
It will be interesting to see the effects of these results on 
astrophysical scenarios. 


IV. SUMMARY 

Low energy (p, 7 ) reactions are studied in a semi- 
microscopic approach in the HE formalism and com¬ 
pared with experiments in the mass region 55-60. Ra¬ 
dial density profiles are obtained using the RME ap¬ 
proach and are folded with the DDM3Y NN interaction 
to obtain semi-microscopic optical potentials. Both the 
real and imaginary depths of the potential are normal¬ 
ized to obtain a good agreement between theory and 
experiment. The S-factors for (p, 7 ) reactions are evalu¬ 
ated in the Gamow window corresponding to 3 GK. We 
have not modified the parameters to fit individual reac¬ 
tions as our aim is to construct a framework for calcula¬ 
tion of astrophysical reactions involving unstable nuclei. 
Rates for important astrophysical reactions calculated 
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FIG. 3: Comparison of astrophysical rates from the present calculation with NON-SMOKER results from Rauscher at al. for the 

reactions (a)®®Ni(p, 7 )®^Cu, (b) ®^Cu(p, 7 )®®Zn and (c) ®®Cu(p, 7 )®‘^Zn. See text for details and the Supplemental Material in Ref. [ 2 HI 
for numerical values of rates. 


in the present approach differ significantly from NON- cutta. 
SMOKER rates. The key feature of our work is that 
we have taken all nuclei in the same footing and same 
methodology has been used for all of them to avoid sys¬ 
tematic error. 
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